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We have identified the optical counterpart of the ULX source P13 in the nearby spiral galaxy NGC 7793. The object is a V 
~ 20.5 mag star, ten times brighter than any other established counterpart of a ULX in nearby galaxies. Medium resolution 
optical spectroscopy carried out in 2008 and 2009 with the ESO-VLT reveals the presence of narrow high order Balmer, 
Hel and Mgll absorption lines indicating a late B type supergiant companion star with mass between 10 and 20 Mq . Stellar 
H/3 and Hell A4686 emission lines are also seen superposed on the photospheric spectrum. We detect different patterns of 
radial velocity variations from the emission and absorption lines over a time interval of one month. The velocity of the 
high order Balmer absorption lines changes by ~ 100 km s -1 while the H/3 and HeIIA4686 emission components vary 
by about the same amount but with a different phasing. Assuming that the observed velocity changes trace the motion of 
the mass-donor star and of the X-ray source implies a mass of the accreting black hole in the range of 3 to 100 Mq with a 
most probable value of ~ 10 to 20 Mq . We expect an orbital period in the range of 20 to 40 days based on the low density 
of the supergiant star. P13 is likely in a short-lived, and thus rare high X-ray luminosity evolutionary state associated with 
the ascension of the donor star onto the supergiant stage. 
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1 The ULX P13 in NGC 7793 

Many galaxies harbour non-nuclear X-ray sources with Lx 
higher than 10 39 and up to several 10 40 ergs _1 , generally 
known as ultraluminous X-ray sources (ULXs). X-ray vari- 
ability shows that the majority are accreting compact ob- 
jects. If the Eddington limit is not violated, or if no strong 
radiation beaming occurs, such isotropic luminosities im- 
ply compact stellar masses up to a few 100 M Q . However, 
presently available stellar evolutionary models at solar me- 
tallicity do not predict compact remnants that massive. Even 
intermediate-mass black holes with M ~ 10 3 M Q have be- 
en suggest ed, based on X-r ay spectral and timing arguments 
(see e.g. Mill er et al. 2004). Determining black hole masses 
is therefore a key issue for understanding the X-ray emis- 
sion mechanism and the evolutionary status of ULXs. 

NGC 7793 is an Sd galaxy member of the Sculptor Gr- 
oup o f galaxies located at a distance of 3.9 Mp c ( Karachentsev 
2003b . ROSAT PSPC observations dRead & Pietschl 1 1 999l) 
have revealed that its X-ray emission is dominated by a sin- 
gle source, P13, located at the southern edge of the galaxy 
and exhibiting evidence of variability over a 5 month time 
interval. Chandra r evisited NGC 7793 in 2003 for 49 ks 
( Pannuti et al] |2006). providing an accurate source position 
at RA = 23:57:51.01, Dec = -32:37:26.62 (J2000) with a 
formal 90% confidence radius of 0.4". A simple power law 
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does not fit well the Chandra ACIS-S spectrum (x 2 =201 
for 181 dof). However, the addition of a soft or of a hot 
blackbody decreases the x 2 to acceptable values (% 2 ~ 185 
for 179 dof in both cases). The corresponding two best fits 
are T = 1.10 ± 0.06, kT = 0.18 ± 0.03 keV and T = 2.14 
±0.28,kT= 1.72±0.13keV. With L x ~ 4 xlO^ergs -1 , 
PI 3 is a bona fide ULX. A blue V-20.5 stellar object falls 
right in the middle of the small Chandra error circle, leav- 
ing no doubt about the identification of PI 3 counterpart, the 
brightest ULX optical counterpart presently known (Pakull 
et al. 2010, in prep). 



2 Optical properties 

In 2008 and 2009 we obtained ESO-VLT observations with 
the a im of deriving the spectral type of the companion star 
eJiftl.Jf measuring the dynamical masses of the binary com- 
ponents. Unfortunately, very few of the scheduled obser- 
vations were eventually executed in 2008 and 2009. Addi- 
tional runs are expected to be carried out in 2010. On 2008 
November 23 and 27, two 40min long spectra were ob- 
tained with FORS 1 and the 600B Grism providing a spectral 
resolution of ~ 600 over the wavelength range 3500-6000A. 
The five 2009 40 min long exposures were obtained with a 
higher resolution of ~ 1600 over the 3700-5200A range us- 
ing FORS2 and grism 1200B. Observations took place on 
November 17 and December 9, 11, 13 and 14. The rectified 
sum of all optical spectra obtained in 2009 shown in Fig.Q] 
reveals many deep unresolved Balmer and Hel absorption 
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Fig. 1 Optical spectrum of P13 obtained by summing all 2009 observations (total exposure time of 3h). The continuum 
has been normalized to unity. Photospheric lines include the high order Balmer absorption lines, some Hel lines and the 
Mgll line. The Hell line is clearly seen in emission. Balmer emission is strong at H/3 and added to the photospheric 
absorption line yields a marked P Cyg profile. 



lines. Strong emission from the Hell A4686 and H/3 lines 
is superposed onto the photospheric absorption spectrum. 
The ratio of the HeIA4471 to Mgl l A4481 lines indicates 
a B8 spectral type (T off - 14.000K: IWalborn & Fitzpatrickl 
1990) while the high optical brightness of the star Mv ~ 
—7 implies a supergiant luminosity class. Our 2008 spectra 
display a strong Balmer jump indicating log g ~ 2 in agree- 
ment with the luminosity class. Masses and radii of such 
stars are no t accurately known. Evolutionary tracks with ro- 
tation from lMevnet & Maederj d2000h predict masses of the 
order of 20 M for our observed Mboi of ~ —8. However, 
masses as lo w as lOMpi have also b een reported for such 
stars (see e.g. iMarkova & Pulsl 120081) . The high optical lu- 
minosity implies a stellar radius of 40 to 60 Rq. Note that 
in spite of its very large optical luminosity, the total emis- 
sion of PI 3 remains dominated by its high energy output. 
The X-ray to stellar bolometric luminosity ratio Lx/Lboi ro 6 
suggests that some X-ray heating effects might be observ- 
able. However, the impact of X-ray heating on the derived 
photospheric T e g is difficult to assess without knowing the 
geometry of the system. 



3 Radial velocities 

We applied the standard ESO gasgano + fors-kit-4.1 pipe- 
line to the 2009 FORS2 data to create wavelength calibrated 
spectra corrected for bias and flat-field. Remaining spectral 
distortions in the near UV region were removed using addi- 
tional arc lines at short wavelengths. The final wavelength 
scale is accurate to 2 kms -1 while the small variation of 
the target position within the slit from one observation to 
the other yielded a scatter of less than 5 km s ~ 1 . 

Fig. Q] shows high order Balmer lines mostly in absorp- 
tion, whereas emission components become conspicuous at 
H7 and below. Using a cross-correlation method we mon- 
itored the velocity of the 5 Balmer lines contained in the 
3750 to 3900A interval. We tested several different tem- 
plates: a) the straight average of all spectra with and without 
velocity correction, b) a B8I spectrum extracted from the 
I Walborn & Fitzpatrickl Jl990h atlas and c) a TLUTSY B- 
STAR06 model dLanz & Hubenvll2007l) with T cff = 15,000 
K and log<? =2.0, degraded to the FORS2 resolution. Re- 
sults are not sensitive to the actual template used or to the 
inclusion of additional absorption lines at slightly longer 
wavelengths such as Ca H&K, He and H<5. Errors on veloc- 
ities were derived from sub-exposures measurements and 
estimated to be of the order of 7 km s _1 . 
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Fig. 2 Variation of the H7 absorption line profile with 
time. Spectra have been shifted in flux for clarity. From top 
to bottom, Nov. 17, Dec. 9, 1 1, 13 and 14. 



In order to monitor the behaviour of Balmer emission 
components, we subtracted from the observed spectra a tem- 
plate photospheric spectrum shifted at the measured veloc- 
ity of the high order Balmer absorption lines. We tested 
again two of the templat es considered for the cross-co rrela- 
tion method, namely the lWalborn & Fitzpatrick (1990) atlas 
spectrum and the TLUTSY model without finding notice- 
able effects on the results. The subtracted spectra only dis- 
play emissions at H7 and H/3. The centroid of the H/3 emis- 
sion line was measured by fitting a Gaussian profile. Un- 
fortunately, the H7 emission is too weak to provide useful 
velocity information. The Hell A4686 line has a broad re- 
solved profile on the top of which a narrower component ap- 
pears on occasion. Since the line does not overlap with any 
strong photospheric feature, we directly measured its po- 
sition on the original spectra by fitting Gaussian profiles to 
the broad or narrow components or using a cross-correlation 
method. Again, all approaches gave quite similar results. 

We find that the velocity of the absorption features chan- 
ges by 105.5 ± 13kms _1 over the 27 day time interval be- 
tween the first and last observation. The velocity shift is also 
visible in the H7 line plotted in Fig. [2] although this partic- 
ular line is clearly contaminated by residual emission. 

Most of the velocity change of the H/3 emission line oc- 
curs between the first and second data points separated by 
22 days (see Fig. [3]) with a total amplitude of 73.4±11.0 
kms -1 . A PCyg profile is patent in several spectra. How- 
ever, the line equivalent width (EW) remains stable at a 
value of ~ —3. 5 A. In contrast, the Hell A4686 line exhibits 
larger profile and velocity variations (see Fig. [4]). Ignoring 
the second observation in which the line was particularly 
broad and weak, the total velocity amplitude of the Hell line 
is 104±15kms~ 1 , i.e. comparable to that of the Balmer ab- 
sorption lines. The line EW varies from —1 to —2.4 A. 

The radial velocity curves plotted in Fig.|5]show a rather 
complex pattern which is hard to interpret considering the 
scarcity of the available data. One can remark, however, that 
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Fig. 3 Variation of the H/3 emission line profile with time. 
A TLUTSY/BSTAR06 model spectrum shifted at the veloc- 
ity of the high order Balmer absorption lines was subtracted 
from the observed spectrum. Spectra are shifted in intensity 
for clarity. From top to bottom, Nov. 17, Dec. 9, 11, 13 and 
14. 
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Fig. 4 Variation of the Hell A4686 emission line profile 
with time. Spectra are shifted in intensity for clarity. From 
top to bottom, Nov. 17, Dec. 9, 1 1, 13 and 14. 



between the single November 17 observation and the mean 
of the December observations, the velocities of the emis- 
sion and absorption components roughly varied in opposite 
directions. Such a behaviour would be expected in a dou- 
ble line spectroscopic binary. The velocity of the absorption 
lines probably reflects to a large extent the orbital motion 
of the mass-donor star. Unfortunately, its smooth variation 
during the 4 day long run in December is not accompanied 
by any correlated change in the H/3 or He II line velocities. 
This hints at an orbital phasing different from 180° between 
absorption and emission lines, or worse, at a mostly ran- 
dom variation. PI 3 might thus be not as clean as one would 
have liked to measure dynamical masses of both compo- 
nents. It is fairly possible that none of the emission lines 
accurately follows the motion of the compact object, nei- 
ther in amplitude nor in phase. Many supergiant stars ex- 
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hibit Balmer emission lines with PCyg profiles produced 
in massive winds. In the case of P13, the difference in ve- 
locity between the peak emission and absorption is ~ 250- 
330 km s~\ fully consistent with the escape velocity from a 
20 M Q star having a 60 R© radius. Therefore, emission from 
an X-ray photoionized stellar wind could well add to the 
light originating from the accretion disc and considerably 
alter the H/3 radial velocity curve. 

Similar considerations apply to the Hell line. Its veloc- 
ity profile might include emissions from the X-ray heated 
companion hemisphere or from various structures of the ac- 
cretion disc. In fact, previous optical spectroscopic obser- 
vations of ULXs cast some doubt on the reliability of the 
Hell A4686 line as a tracer of t he motion of the accreting 
component ( Roberts et al. 2010h . Most ULXs are likely to 
undergo noticeable X-ray heating effects with an irradiated 
to unirradiated optical flux ratio closer to that seen in galac- 
tic low-mass X-ray binaries than shown by classical high- 
mass X-ray binaries (HMXRBs) dCopperwheat et al.f 2005: 
Grise et alj|2008l) . This could imply, for instance, that the 



Hell A4686 line does not follow the velocity of the accret- 
ing component but rather that of the disc bulge at the point 



where the stream of matter escaping L I impacts the edge of 

_ ]|2006). Moreover, although 



the disc (see e.g. iPearson et al 



for the first time in a ULX we detect photospheric absorp- 
tion lines, in other ULXs spectral sign atures from the com - 
panion stars have not yet been found ([Roberts et al . 2010). 

However, in the case of P13, the roughly opposite ve- 
locity variations of the Hell + H/3 and high order Balmer 
absorption lines suggest that these lines may still be used to 
constrain to some extent the dynamical parameters. In addi- 
tion, the moderately high Lx/Lboi of ~ 6 is similar to that of 
the X-ray luminous HMXRB LMC X-4 in which the Hell 
A4686 line does trace the motion of the neutron star (Kx ~ 
KM<,TT: lKellev et al.lll983h . 



Adopting the Roche lobe overflow formalism of EggletonJ 



(119831) yields orbital periods of about 20 to 40 days for a 
large range of mass ratios and for stellar radii between 40 
and 60 R©. The fact that the total velocity amplitudes of 
Balmer absorption and Hell lines are both larger than ~ 
100 km s -1 implies a mass ratio q = M opt /Mx in the range 
of 0.2 to 3.2. Values of Kx ~ K opt ~ 50 km s -1 are reached 
for i ~ 30°. A mass ratio close to 1 is supported by the 
nearly equal amplitude of the absorption and emission lines 
radial velocity curves and therefore hints at a black hole 
mass of ~ 10 - 20 M Q . Taken at face value, our observa- 
tions do not support the presence of an intermediate mass 
black hole in PI 3. 



4 Discussion and Conclusions 

Determining whether ULXs harbour stellar mass or in some 
cases intermediate mass black holes can only be settled by 
measuring orbital periods and dynamical masses. However, 
the optical faintness of most ULXs is a challenge for me- 
dium resolution spectroscopic monitoring, even using 8m- 
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Fig. 5 Radial velocity curves of the high order Balmer 
absorption lines (black filled squares) of the H(3 emission 
(open squares) and of the He A4686 emission line (red 
filled squares). The horizontal line shows the velocity of the 
nearby interstellar gas. 



class telescopes. P13 is the first ULX in which optical spec- 
tral signatures of the mass-donor star are detected. These 
lines provide both a reliable indication of the spectral type 
and luminosity class of the companion star and at the same 
time means to measure accurate orbital velocities. The high 
optical luminosity and the late B spectral type of the com- 
panion star clearly demonstrate that the star has left the main 
sequence and is now rapidly evolving towards a red super- 
giant. This evolutionary stage c orresponds to the most X -ray 
luminous phases of any ULX dRappaport et al.l 120051) and 
should last no mor e than ~ 10 6 years and m ay be as short 
as a few 10 4 years ( Patruno & Zampieril 2008 ). In compar- 
ison, the life time of ULXs in which the mass-donor star is 
a sub-giant is much longer (~ 10 7 years). Population mod- 
els indeed confirm that systems like PI 3 are rare and shoul d 
account for only ~ 5% of all ULXs (Rappaport et al ]l2005h . 
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